Phonon Rabi-assisted tunneling in diatomic molecules 
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We study electronic transport in diatomic molecules connected to metallic contacts in the regime 
where both electron-electron and electron- phonon interactions are important. We find that the 
competition between these interactions results in unique resonant conditions for interlevel transitions 
and polaron formation: the Coulomb repulsion requires additional energy when electrons attempt 
phonon-assisted interlevel jumps between fully or partially occupied levels. We apply the equations of 
motion approach to calculate the electronic Green's functions. The density of states and conductance 
through the system are shown to exhibit interesting Rabi-like splitting of Coulomb blockade peaks 
and strong temperature dependence under the interacting resonant conditions. 
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A significant current effort in nanoscopic systems is 
the study of electron transport in natural and quantum- 
dot molecules. Much of the interest lies in being able 
to investigate different regimes of competing electron- 
electron and electron-phonon interactions. It is typically 
the case, due to the spatial confinement, that electron- 
electron interactions (EEI) play a more important role 
than electron-phonon interactions (EPI) in determining 
electronic transport properties in low-dimensional sys- 
tems. 

Different geometries of quantum-dot molecules (con- 
structed with interconnected quantum dots) have been 
studied in the literature, and the role of phonons on 
electron transport has been analyzed in these systems^ 
For instance, it is known that phonons are a relatively 
weak perturbation, responsible for the broadening of 
Coulomb blockade peaks in the conductance and for the 
appearance of satellite features in the nonlinear transport 
regime. 1 

More recently, the field of "molecular electronics,"^ 
where electrons and/or holes are injected directly into 
molecules attached to metal electrodes, has seen intense 
activity and progress4£ It is interesting to note that EPI 
become more important in molecular electronics, since lo- 
cal molecule deformations produce significant electronic 
level shifts, as has been observed in experiments.- In fact, 
vibrational and torsional modes play prominent roles in 
electron transport, producing sidebands in the voltage- 
dependent differential conductanc e) 7 ' 8 ' 9 and/or polaronic 
shifts of the electronic levels. 10 Furthermore, in a molec- 
ular system with discrete electronic energy levels, vibra- 
tions produce important effects when the energy of the 
vibrational modes matches the energy difference between 
electronic levelsiii As a result, EPI provides relaxation 
mechanisms (inelastic scattering) that affect the conduc- 
tance of the system^ 

The simplest model of EPI is perhaps the indepen- 
dent boson models where localized electrons interact 
with a phonon system. Phonons introduce a shift of elec- 
tronic levels and create a series of phonon replica peaks 



in the density of states ( DOS) of the electronic sys- 
tem. In a double barrier heterostructure interacting with 
phonons, as electrons have access to an energy continuum 
in the region outside the barriers, the inelastic scattering 
strongly affects the resonant tunneling regime through 
the heterostructureii^ The effects of EPI on a double- 
level quantum dot have also been studied recently^ al- 
though the combined effects of inelastic scattering and 
EEI were not considered. 

In order to study a more realistic molecular system it 
is important to include both interactions — EEI and EPI 
- simultaneously. The interplay between the compet- 
ing interactions is likely to result in unique conditions for 
phonon emission and absorption, as well as in unexpected 
polaron behavioriiS In this work we study the effect that 
EEI and EPI have on charge transport through a di- 
atomic molecule, envisioned as two atomic sites (or quan- 
tum dots) directly coupled to leads, as shown schemati- 
cally in Fig.Q] We study this system using the equations 
of motion method, which allows us to obtain the DOS and 
electronic occupation as well as the conductance through 
the system. We exploit the fact that the EPI strength 
A is small compared to the phonon energy, and thus in- 
clude self-energy terms up to second order in this param- 
eter. Standard considerations, similar to those applied 
in the Hubbard approximation^^ 7 " are used to evaluate 
the equations. A significant result is the identification 
of occupation-dependent resonant conditions for phonon 
absorption and emission in the presence of Coulomb re- 
pulsion. More importantly, we find a unique type of Rabi 
splitting in the DOS from the mixing between a doubly 
occupied low-energy level (boosted by the Coulomb re- 
pulsion) and a higher-energy state. This Rabi splitting 
is mediated by thermally regulated phonon emission and 
absorption in the molecule. The effect is shown to dra- 
matically modify the transport properties of the system 
since the resulting polaron formation competes with reso- 
nant tunneling. The effect is made remarkably noticeable 
even for weak EPI since the phonon-assisted transport is 
magnified by the virtual emission and absorption pro- 
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FIG. 1: Schematic representation of the model system. The 
electron-phonon interaction connects the local sites (e a < eg) 
via phonons of frequency u>o with coupling constant A. 



cesses in the interacting resonant regime. 

For concreteness we consider a two-level diatomic 
molecule with local energies e° and e|g (we assume e° > 
e°), as shown in Fig. ^ Each dot or atomic site is 
connected independently to two external current leads. 
The system as described can be mapped to that stud- 
ied in Ref. 0| (there for the one-electron case) and is 
designed to model various experimental geometries. The 
total Hamiltonian is written as Ht — H m ol + Hi ea d s + 
H m oi-ieads- Each lead is modeled as a semi- infinite tight- 
binding chain, Hi ead = J2a.<j',j> tc ]'a c j^' where the 
site index sum is over nearest neighbors, ct (cj a ) cre- 
ates (annihilates) a fermion at the j-th site with spin 
a, and the left (right) lead is defined for < — 1 

(J'jj > !)• The Hamiltonian for the molecule is given 
by H mo i = H e i + H p h + Hei^ph, 
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where b^ (b) creates (annihilates) a phonon with energy 
fkdo, ei = e" — eV g , and i = a,(3. The gate voltage V g 
controls the particle number by shifting both localized 
energies with respect to the Fermi energies of the left 
and right electrodes. We assume here the same gate 
on both sites, as is likely the case in molecules^ Fi- 
nally, H mo i- leads = X^.™.- 3 Uc\ a Cj a + H.c, connects 

3=1,-1 

the diatomic molecule to the leads, where we will con- 
sider t a ,tp <C t. Away from the Kondo regime the effect 
of the leads is to broaden the energy levels of the dots 
through the tunnel couplings t a and tp. (Ref. 0) 

To determine the dynamics of the electrons in the 
molecule, we calculate the local retarded Green's func- 
tion, together with the nonlocal functions, 
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The latter are needed since they describe electron prop- 
agation due to the EPI, and are associated with the ab- 
sorption (G^) and emission (Gp a ) of phonons. The cor- 
responding equations of motion up to second order in the 
EPI coupling strength A, yield (i = a,0), 
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Gp a (e) = 
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where the expressions for the corresponding self-energies 
are given by 
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e — e Q — U{n a ) — tkoo + iT(e — Hluq) 

We have considered here the paramagnetic case with 
(nia) = {fiia) = {rii). The broadening due to the leads is 



given by Ti(e) — 2ntfp(e), where p(e) is the DOS in the 
leadsi We take T a = Tp = T(e), for simplicity^ 

The nonlocal Green's functions, Eqs. |[BJ and Q, have 
simple physical interpretations. At low temperature T, 
{bib) w thus Gp a -» as {np) — > 0, i.e., there is no 
phonon emission if the (3 level is empty. It is also inter- 
esting to note that if the dot a is not completely full, 
{n a ) < 1, the process of phonon absorption described by 
G a p is possible (even at small T). 

In the limit U = the local Green's functions Gu 
have two poles, one at e$ and another at the pole of the 
self-energy E,. In this case there is a single phonon reso- 
nance condition for both levels (i = a, /?), achieved when 
the phonon energy tkoo matches the energy difference be- 
tween the two localized electron energies Ae = ep — e a , 
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FIG. 2: (Color online) Total density of states (DOS) vs en- 
ergy and gate voltage for a two-site molecule as in Fig. Q 
In units of Ae, the parameters are: A = 0.2, fcsT = 0.015, 
t a — tp — 0.1, U — 0.4, and hwo = 0.6. The Fermi level is set 



at zero energy; e a — 1 and ep = 2 for V g = 0. For V g ~ 
2.4, the DOS exhibits a clearly split feature at 
produced by the strong phonon-assisted level coupling. 
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inset shows the DOS at V g = 2, indicated by the dashed line 



in the main figure. 



i.e., at Ae = Hluq. Thus, the resonance effectively couples 
both sites. Such a phonon resonance condition on trans- 
port has been recently explored by Tasai and Eto^ They 
find a sharp dip in the conductance as the result of de- 
structive interference between bonding and antibonding 
states (see the dashed line in Fig. 2). We will see below 
that electron repulsion greatly affects this behavior. 

In the limit A = we recover the Hubbard I 
approximation^ As expected, each local Green's func- 
tion has poles at and e; + U. In contrast, when U =/= 0, 
we find different resonant conditions for a and j3. The 
Coulomb repulsion requires extra energy for electrons to 
tunnel into a fully or partially occupied state. This ex- 
tra energy depends on the occupation fraction of the two 
sites, which can be controlled by the gate voltage. The 
self-consistent charge of each site is obtained by integrat- 
ing the DOS, pa — (— l/7r)lmGji(e), for different gate 
voltages. With appropriate parameter values we find that 
the effects brought about by the EPI and EEI are em- 
phasized under resonant conditions. Hereafter all energy 
quantities are given in units of the level spacing Ae. 

Figure shows the total (for both sites) DOS for the 
molecular system as a color map for different values of 
gate voltage V g and energy w, measured with respect to 
the lowest level e Q . The temperature is set at fc^T = 
0.015, U = 0.4, and hw = 0.6. [The DOS is normalized 
to obtain the correct charge number for each dot.] For 
low-voltage values, the molecule is empty and the DOS 
shows two main features at the noninteracting energies 
uj ~ e a , e@, slightly shifted and broadened by the coupling 
to the phonons and leads. As the gate voltage increases, 
the lowest level approaches the Fermi level, and as its 
occupancy grows a unique feature in the DOS develops. 
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FIG. 3: (Color online) Conductance vs gate voltage for 
fcsT/Ae = 0.015, 0.10, and 0.15. Temperature has a strong 
effect on the resonance condition, enhancing the Rabi splitting 
of the third Coulomb blockade peak (V g — e Q ~ 1) at higher 
T and exposing the shoulder in second peak. The inset shows 
inelastic contributions for the conductance. The highest peak 
is due to phonon emission. The dashed line shows results for 
U — and wo = 1. Curves offset for clarity. 



This is shown in Fig. [21 at energy e a + U » 0.4 for V g > 
1. Notice that at the same time, a weak phonon-related 
feature (a "phonon replica") appears at energy ~ e Q + 
Hujq = 0.6, as can be seen from the expression for the 
self-energy, Eq. @ . [This feature persists weakly even at 
low V g at ~ 0.5, due to the finite A and T values in Fig. 

F° r V g > 1.4, the lowest level is almost fully occupied, 
(n a ) ~ 1, and the resonant absorption/emission condi- 
tion in Eq. @ is achieved when e a + U{n a ) + Hujq ~ ep. 
The increasing (n a ) occupation shifts the phonon replica 
feature, moving it closer to ep. The resonant condition 
results in the effective mixing of the e Q + U and ep levels 
due to phonon-assisted transitions. The mixing produces 
a near degeneracy and, as a consequence, an effective 
phonon Rabi splitting of spectral features appears. This 
is the origin of the doublet appearing at u> — e a ~ 1 with 
nearly equal-size peaks in the range V g ~ 1.6-2.4. Note 
that the lower level participates in the process as a re- 
sult of the Coulomb repulsion between electrons. Notice 
further that the resonant condition disappears once the 
higher level becomes fully (doubly) occupied (V g > 2.4 in 
Fig. 13). The DOS returns then to the standard Hubbard 
peaks at e a + U and ep + U dominating the spectrum of 
the molecule (weak phonon-related features near and 
0.8 are present for finite T and A). As we will show be- 
low, the Rabi resonance in the presence of EEI affects 
also the conductance through the system, providing an 
experimentally accessible signature of the effect. 

Figure shows the conductance G vs gate voltage 
for different temperatures (G oc | 2 where G\x oc 
t 2 a G aa + t 2 p Gpp + t a [G a p + Gp a ]tp). As the gate volt- 
age increases, the conductance exhibits the anticipated 
Coulomb blockade (CB) peaks. Notice that the first two 
are associated here with e Q , and show nearly full e 2 /h 



4 




0.96 0.98 1.00 1.02 1.04 



FIG. 4: (Color online) Conductance through the molecule 
at around the third Coulomb blockade peak as a function 
of interaction U and gate voltage V g for ksT = 0.015Ae. 
The inset indicates phonon emission and absorption processes 
involved. 

conductance (limited by finite temperatures). The dom- 
inant effect of higher temperatures is to produce impor- 
tant changes in the dot occupancy due to phonon emis- 
sion (occupancy changes due to thermal excitation of 
electrons are negligible). The third CB peak appearing 
at V g — e° a ~ 1 is clearly split, indicating the phonon- 
mediated transitions between the a and f3 levels. The 
temperature dependence of this peak is pronounced, as 
the thermal variation in electron and phonon occupations 
will not only make it weaker but will also rapidly detune 
the resonance condition. A higher phonon presence re- 
sults in stronger EPI and a drop in the conductance, a 
typical signature of interference between resonant tunnel- 
ing and polaron formation processes. The splitting of the 
third CB peak increases linearly with the EPI strength 
A, as expected in Rabi splitting phenomena. Notice that 
once the j3 level is nearly full, the effect disappears, and 
the Coulomb blockade peak at V g — e° ~ U + Ae = 1.4 is 



essentially T independent. The appearance of a shoulder 
in the second CB peak for higher temperatures is also a 
consequence of EPI, but one that is inelastic in nature: 
The inset of Fig. [3] shows the contribution of inelastic 
processes (oc G a p + Gp a ) to the conductance. The high- 
est peak at V g — e° a ~ 0.4 is due to phonon emission 
processes: electrons thermally excited to level tp, (above 
the Fermi level here) fall to level e a + U by the emission 
of phonons (the inset of Fig. ^illustrates levels involved). 
These processes do not enhance the conductance of the 
resonant level, but rather reduce it, as one sees in the 
main panel. The suppression is produced by the destruc- 
tive interference of the two different conducting processes 
that electrons undergo. 

The resonant condition of phonon-assisted tunneling 
through the molecule in the presence of EEI is a func- 
tion of the interaction parameter U and the occupation 
of the molecule. In order to explore the dependence of 
the Rabi splitting on interaction parameters we show in 
Fig. 0] a color map of the conductance as function of 
U/Ae and (V g — e°)/Ae. It is clear that the conductance 
returns to the non-phonon-assisted resonant single CB 
peak whenever U/Ae is far below or above the resonant 
value U ~ Ae — huiQ — 0.4. The lowest level is nearly dou- 
bly occupied in this V g regime, i.e., (V g — e° )/Ae ~ 1. At 
the resonance, U ~ 0.4, the conductance Rabisplits into 
two peaks that have the same height and demonstrates 
the effect of polaron formation on the conductance. 

We have shown that the competition of EEI and EPI in 
a diatomic molecule produces unexpected Rabi splitting 
phenomena in the DOS with observable effects on the 
conductance. This phenomenon involves states produced 
by the Coulomb repulsion, and it is enhanced at higher 
temperatures, a direct consequence of the thermal nature 
of the phonon bath involved. 
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